Exact calculation of thermodynamical quantities of the integrable t — J model 
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The specific heat and the compressibility for the integrable 
t — J model are calculated showing Luttinger liquid behav- 
ior for low temperatures. A Trotter-Suzuki mapping and the 
quantum transfer matrix approach are utilized. Using an alge- 
braic Bethe ansatz this method permits the exact calculation 
of the free energy and related quantities. A set of just two 
non-linear integral equations determining these quantities is 
studied for various particle densities and temperatures. The 
structure of the specific heat is discussed in terms of the ele- 
mentary charge as well as spin excitations. 
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The t — J model has been studied intensively in re- 
cent years because of its importance as one of the most 
fundamental systems of strongly correlated electrons and 
its possible relevance for purely electronic mechanisms for 
high-T c superconductivity . The model describes the 
nearest-neighbor hopping of electrons with spin-exchange 
interaction. The effect of a strong repulsive on-site 
Coulomb interaction is modeled by the restriction of the 
Hilbert space to states without doubly occupied lattice 
sites. The one-dimensional Hamiltonian reads 

+ J ^2( S J S j+T- - %%+l/ 4 ) (!) 

i 

where the projector V = Ylj(^ ~ n j] n jl) ensures that 
double occupancies of sites are forbidden. At the super- 
symmetric point 2t = J the system was shown to be 
integrable by the well-known Bethe ansatz. The 
ground state and excitation spectrum were investigated 
by Bares et al. ||. Critical exponents were calculated by 
finite-size scaling and conformal field theory dj7|. 

According to the pioneering work of Takahashi M ther- 
modynamical properties of general integrable systems are 
described by an infinite set of coupled integral equations 
reflecting the existence of infinitely many different rapid- 
ity patterns. For the t — J model such integral equations 
were formulated in PJlCfl, but not solved. We like to 
point out two shortcomings of the traditional approach. 
First, the method consists of a rather direct evaluation of 
the partition function by taking into account all excited 
states of the Hamiltonian. The excitations are derived 
on the basis of the so-called string conjecture which de- 



scribes the solutions of the Bethe ansatz equations. How- 
ever, the validity of the string conjecture is quite con- 
troversial. Second, the calculation of the asymptotics of 
correlation functions is not possible within the traditional 
approach. 

Our approach takes a different route following pi 12 
to overcome the mentioned problems. We do not directly 
study the Hamiltonian at finite temperatures. Instead, 
we employ a convenient Trotter-Suzuki mapping ||i"lf to 
an exactly solvable two-dimensional classical model. Fur- 
thermore, the corresponding thermodynamical quantities 
are expressed by eigenvalues of the so-called quantum 
transfer matrix. The largest eigenvalue, for example, 
directly yields the free energy. The next-largest eigen- 
values provide the correlation lengths. This represents 
the important advantage compared with the traditional 
thermodynamical Bethe ansatz requiring all eigenvalues 
of the Hamiltonian. The new method using the Trotter- 
Suzuki mapping as well as the quantum transfer matrix 
approach has been applied to several quantum systems 
jfl] {L6| notably the Hubbard model in Q . 

The t — J quantum chain belongs to a family of inte- 
grable fermionic systems which are related to the classi- 
cal Perk-Schulz model a multi-component general- 
ization of the well-known six-vertex model. The detailed 
derivation of the quantum transfer matrix and the study 
of the corresponding eigenvalue problem by the quantum 
inverse scattering method will be described in a separate 
publication jl1||2(J. Thus, we confine ourselves to the 
main points. 

The Perk-Schultz model is exactly solvable as the 
Boltzmann weights satisfy the Yang-Baxter equation. 
This in turn implies the commutation of all row-to-row 
transfer matrices for arbitrary spectral parameters u, v: 
T(u)T(v) = T(v)T(u). Consequently, the Hamiltonian is 
integrable. We have 

H = -^-lnrfuY =Q with T(u) =T R e uH+a{u2 \ (2) 

where Tr is a right-shift operator. By means of the sub- 
stitution u — —fi/N (/3 being the inverse temperature) 
we find 



(T(-p/N)T(-p/N)) N/2 = e -PH+o[i/N)^ 



(3) 



where T is the row-to-row transfer matrix of the Perk- 
Schultz model after a 90° rotation. It satisfies (||) with 
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T, Tr replaced by T, T R . The partition function of the 
quantum system 



Z = lim Trace(T(-/3/N)T(-0/N)) 



N/2 



(4) 



is identical to the partition function of the inhomoge- 
neous Perk-Schulz model with alternating rows fl9fl . The 
technically more convenient column-to-column transfer 
matrix of such a system is often referred to as the quan- 
tum transfer matrix (QTM). This matrix is shown to 
be embedded into a commuting family 19 2(1 of stag- 



gered transfer matrices given by alternating products of 
unrotated and rotated vertices. We still have the same 
intertwiner as in the case of the homogeneous lattice. In 
contrast to the work in |l7j for the Hubbard model the 
integrability is preserved within the present approach. 
Thus, the eigenvalues of the QTM are analytic functions 
of the spectral parameter. This fact represents the im- 
portant advantage leading to simple calculations which 
can be generalized to other models. The free energy can 
be calculated by the largest eigenvalue A of the QTM 
/ = — (logA)//3 and the next-leading eigenvalues yield 
the correlation lengths. 

We treat the quantum transfer matrix by using an al- 
gebraic Bethe ansatz JlSj ] . The eigenvalues A(w) of T(v) 
are shown to be analytic functions of the spectral param- 
eter v with explicit representation 



A(v) = X 1 (v)+X 2 (v) + X 3 {v), 



(5) 



where 
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V — Vk 



(v — iu — i)(v + iu) 



(v + iu + i){v — iu) 



N/2 



N/2 



(if 



v - Wj +i ^ 



and Vj and Wj are Bethe ansatz rapidities, /1 denotes 
the chemical potential measured from the lower edge of 
the band. The defining relations for the rapidities, the 
so-called Bethe ansatz equations are equivalent to the 
analyticity of A(u), i.e. the absence of poles. We use 
this analyticity to determine the largest eigenvalue by a 
finite set of non-linear integral equations [^9|,^0|. This 
also allows for taking the limit N — + 00 analytically. We 
define the following auxiliary functions (with arbitrary 
but fixed < 7 < 1) 



b(x) 
c(x) 



Ai(x + 17) 



X 2 (x + 17) + A 3 (a; + 17) ' 

Xi(x)X 2 (x) 

A 3 (a;)(Ai(a;) + A 2 (a;) + A 3 (x)) 



Using this definition and the analyticity property of the 
eigenvalue A(i>) a set of functional equations for b and c is 
derived. This set can be transformed into integral equa- 
tions by standard applications of the Fourier transform. 
Such an approach should be applicable to all integrable 
models. After some lengthy calculations (20| we obtain 
the following relations for the two auxiliary functions 

logfr(x) = -2TTf3V b (x + ij) + /3fi - 

* 6 *log(l + 6)| *6 * log(l + c)\ x+il , 

\ogc{x) = -2tt/3^ c (x) + 2(3(1 - 

25R(* b * log(l + 6)U +i7 ) - * c * log(l + c)U, 

with 

$ b (x) = l/(2wx{x - »)), ^c(x) - 2/(2^(x 2 + 1)), 

where * denotes the convolution f*g = J f(x — y)g{y)dy 
taken at the indicated arguments x, x + -17 and x + 2i^. 
The solution of the integral equations provides the largest 
eigenvalue via: 



log A = -logc(O) + 2/3/1. 



(6) 



In order to calculate derivatives of the thermodynami- 
cal potential one can avoid numerical differentiations by 
utilizing similar integral equations p(| guaranteeing the 
same numerical accuracy as for of the free energy. 

In the low-temperature and low-density limit (small 
fi and T/fi) we obtain an analytical expression for the 
grand canonical potential / = — (log A) / (3 



dx 

7T/3 3 / 2 



log (1 + e 



(7) 



in agreement with the numerical results. Asymptotically 
we find 

/ = ~ M 3/2 " ^= T 2 + o{^\ T 2 /VM), (8) 

implying a particle density n = —df/d/i and entropy 
S = -df/dT: 



2 2 
n = -^/JI, S=—T. 
it 6n 



(9) 



Moreover, due to predictions of Luttinger liquid the- 
ory and conformal field theory we expect for the low- 
temperature asymptotics 



(10) 



where v StC and c 



are the velocities and central charges 
(c^c = 1) for the elementary spin and charge excitations. 
For small particle densities we have u SjC = im ||, thus 
(||) and (|h]) are consistent. 
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On the other hand, the simple high-temperature limit 
for finite \i immediately leads to 5 = ^3 (and n = 2/3) 
as is expected by counting the degrees of freedom per lat- 
tice site. Therefore, the analytical benchmarks for T — > 
and T — > oo imply the consistency of our approach. 

Now we consider various thcrmodynamical quantities 
at intermediate temperatures by solving the integral 
equations numerically. Figure [l] shows the specific heat as 
a function of T for different fixed particle densities. First 
of all, we note a linear temperature dependence at low T. 
According to conformal field theory (|l^) the coefficient is 
given by ir(l/v s + l/w c )/3. Our numerical data are con- 
sistent with this expression. Furthermore, we observe two 
maxima with changing dominance for increasing particle 
density n. The nature of this structure can be under- 
stood from the elementary excitations of the system. In 
the groundstate the particles are bound in singlet pairs 
with binding energies varying from to some density de- 
pendent value. There are two types of excitations. First, 
there are charge excitations due to energy-momentum 
transfer onto individual pairs. Second, there are excita- 
tions due to the breaking of pairs. The latter excitation 
is of spin type at lower excitation energies, but changes 
the character at higher (density dependent) energies to 
charge type as it describes the motion of single particles. 
Therefore, the first and second maximum at lower densi- 
ties (figure |l].a) are caused by charge excitations due to 
pairs and single particles, respectively. At higher densi- 
ties (figure |l].b) the maximum at lower temperatures is 
dominated by excitations of pairs whereas the second one 
at higher temperatures is caused by spin excitations. For 
increasing concentration the spin contribution becomes 
dominant as the charge excitations freeze out. This is in 
accordance with the limiting case n = 1 leading to the 
spin- 1/2 Heisenberg chain. The missing spin structure 
in the specific heat at low and intermediate densities is 
found at quite low temperatures shown in figure |[ 

It is worthwhile to compare these results with the find- 
ings for the Hubbard model investigated by the tradi- 
tional thermodynamical Bethe ansatz The struc- 
ture found in the specific heat is explained by spin and 
charge excitations which do not change their character 
in contrast to the t — J model. For certain densities a 
low-temperature charge peak was found which is caused 
by single particle excitations. A charge peak at higher 
temperatures appears because of excitations due to dou- 
bly occupied lattice sites with energies of the order U for 
large Coulomb interaction. 

To conclude our investigation we present numerical re- 
sults for other thermodynamical quantities. Figure || 
presents the entropy and the compressibility K — dn/d/i 
for various particle densities. Note the divergent low- 
temperature compressibility for particle densities n — > 
and n — > 1. 

In summary, we have derived eigenvalue equations for 
the quantum transfer matrix of the t — J model. This 
approach permits the exact analytical as well as numer- 
ical calculation of thcrmodynamical quantities. Instead 



of solving an infinite set of integral equations - as is nec- 
essary in the traditional thermodynamical Bethe ansatz 
- we have to solve integral equations for only two func- 
tions. We have considered analytically certain low and 
high-temperature limits verifying the values predicted by 
conformal field theory. Moreover, the case of intermedi- 
ate temperatures was treated numerically. As shown, the 
specific heat and compressibility display an interesting 
behavior in dependence on particle density and temper- 
ature. 
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